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ABSTRACT 

We analyze the ionization state of the magnetohydrodynamically turbulent 
protoplanetary disks and propose a new mechanism of sustaining ionization. 
First, we show that in the quasi-steady state of turbulence driven by magnetoro- 
tational instability in a typical protoplanetary disk with dust grains, the amount 
of energy dissipation should be sufficient for providing the ionization energy that 
is required for activating magnetorotational instability. Second, we show that in 
the disk with dust grains the energetic electrons that compose electric currents 
in weakly ionized gas can provide collisional ionization, depending on the actual 
saturation state of magnetorotational turbulence. On the other hand, we show 
that in the protoplanetary disks with the reduced effect of dust grains, the tur- 
bulent motion can homogenize the ionization degree, leading to the activation of 
magnetorotational instability even in the absence of other ionization processes. 
The results in this Letter indicate that most of the regions in protoplanetary 
disks remain magnetically active, and we thus require a change in the theoretical 
modeling of planet formation. 

Subject headings: protoplanetary disks - - planetary systems: formation - 
instabilities — MHD — turbulence — accretion, accretion disks 



1. Introduction 



Protoplanetary disks are supposed to be the sites of planet formation. In the standard 
models of planet formation, the effect of the magnetic field is not seriously considered be- 
cause the ionization degree is supposed to be very small at radii where the planet formation 
should take place. The low ionization is due to the high density of the midplane of the disk, 
where the recombination rate is too high to maintain ionization only through cosmic rays. 
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However, if the ionization degree is sufficiently high, the magnetorotational instability (MRI) 
is important and the disk becomes turbulent (Balbus & Hawley 1991, 1998; Stone et al. 
2000). The combination of Maxwell and Reynolds stresses in magnetohydrodynamic turbu- 
lence provides an efficient mechanism for the radial transport of the angular momentum in 
the disk, which results in gas accretion onto the central star through the disk. Therefore, 
the determination of ionization degree in realistic protoplanetary disks is very important and 
remains the subject of extensive research (Sano et al. 2000; Glassgold et al. 2000; Fromang 
et al. 2002; Salmeron & Wardle 2003; Matsumura & Pudritz 2003). The existence of dust 
grains in the disk has a drastic effect on the ionization degree, because the recombination 
on the dust grains is very efficient in the high-density regime. Most of the studies on the 
ionization degree with the effect of dust grains show that there is a dead zone in which gas 
effectively decouples with the magnetic field. The location of the dead zone spans from a 
few tenths of AU to tens of AU from the central star depending on the (column) density 
distribution of the disk. 

The existence of the dead zone should have important implications for the evolution of 
protoplanetary disks and the formation of planets (e.g., Gammie 1996, Matsumura & Pudritz 
2005) and their migration in the disk (e.g., Terquem 2003, Laughlin, Steinacker, & Adams 
2004, Nelson &Papaloizou 2004), which requires more detailed analysis on the validity of its 
existence and its precise location. In this Letter, we reanalyze the ionization degree taking 
into account additional sources of ionization that are not considered previously. 



Ionization Degree The linear stability analyses and nonlinear simulations show that the 
MHD turbulence in the disk is driven by the MRI if the magnetic Reynolds number is greater 
than about unity: Reyi = v\/(Qrj) > 1 (Sano, Inutsuka, & Miyama 1998; Fleming, Stone, 
& Hawley 2000; Sano & Stone 2002a,b), where va is the Alfven velocity and the resistivity 
7] can be expressed as a function of the electron number density x e and temperature T: 
rj ks 2 ■ 10 2 T 1 / 2 /x e . Thus, the condition is expressed as 



where C s denotes the sound speed, H is the scale height of the disk, and (5 = 2Cl/v\ is the 
so called plasma (3. In order to satisfy Reu > 1, the number fraction of electrons should be 
larger than x e = 10~ 13 at r « 1AU for typical protoplanetary disks with (3 = 10. Numerical 
simulations on the saturation level of the nonlinear development of MRI have shown that a 
typical value of (3 ranges from 10 to 100 (e.g., Hawley, Gammie, & Balbus 1995, Sano et al. 



2. Ionization Required for MRI 
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2004). In this Letter we adopt (3 = 10 as a fiducial value that corresponds to B pa 6 Gauss 
for n = 10 15 cm- 3 and T « 3 • 10 2 K. 

Scaling Relation in High Density Regime with Dust grains The fractional number 
densities of various charged species in the protoplanetary disks with dust grains are essentially 
determined by a single parameter that is the ionization rate ( divided by the number density 
n pa 2n H2 and calculated by Sano et al (2000). Figure 1 plots the fractional number densities 
of three important species as a function of (/n. The range of the fractional number density 
of electrons required for MRI, x e ~ 10~ 13 , corresponds to the normalized ionization rate 
C/n pa 3 • 10~ 28 sec -1 cm 3 , i.e., 

C « 3 • 10- 13 ( " - ) sec" 1 . (2) 
V10 15 cm- 3 / v ' 

This shows that for n = 10 15 cm -3 MRI requires an ionization rate 30,000 times higher than 
the standard cosmic-ray rate (Ccr ~ 10~ 17 , see, e.g., Spitzer 1978). In the following, we 
consider whether this ionizing source is available or not. 

3. Feedback from MHD turbulence 

The column density of typical protoplanetary disks is so high that the external ionizing 
radiation (UV and X-rays) cannot penetrate deep into the midplane. We should consider the 
other mechanisms of ionization. Sano & Inutsuka (2001) have shown that the energy dissi- 
pation rate Q[erg cm -3 sec -1 ] in the turbulence driven by MRI in the disk can be calculated 
by the time and spatial average of the r0-component of the stress tensor: 

Q = t ((^ v * - ^r)) w °- 03 n « B2 »> ^ 

where the double angle brackets denote time and spatial average (see also Table 4 of Sano et 
al. 2004). Most of this energy dissipation is supposed to be used for heating of the gas. In a 
quasi-steady state condition, the input thermal energy is eventually converted to radiation 
that escapes from the disk. Here we propose that some fraction of the energy dissipation 
should be used for ionization. 

Energy Budget for Sustaining Ionization First we consider whether the energy input 
available in the MHD turbulence can provide sufficient energy for the ionization required for 
MRI. The ratio /ionize of the energy required for the ionization to the energy available in the 
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MRI-driven turbulence can be calculated as 

^ionize \ / , \ 

/ionize = — ^— = U.Ud ^ 1015cm _ 3 ; 1^— ^~ J ^ -gj {Y3SeVJ ■ [ ] 

This is sufficiently small for typical conditions in protoplanetary disks even at 1 AU. There- 
fore, if the energy dissipation in MHD turbulence can be used for ionization of the gas, the 
positive feedback loop of MRI can maintain MHD turbulence in the protoplanetary disks as 
schematically shown in Figure 2. 

4. Available Ionization Processes 

The next question is about the microscopic processes of ionization: What is the process 
responsible for ionization in weakly ionized dusty gas? 

Thermal Ionization Many numerical studies have shown that the saturated state of 
MRI-driven turbulence is the so-called high-/? plasma: (5 = 2C 2 /t> 2 v > 10; i.e., the magnetic 
energy is smaller than the thermal energy of the gas. Therefore, even a rapid conversion of 
all the magnetic energy into thermal energy (via, e.g., magnetic reconnection) can raise the 
temperature of the gas only slightly. Thus, the magnetic dissipation does not seem to result 
in thermal ionization. 

Note, however, that this argument requires great care because the thermal ionization 
of some alkali metals takes place at quite low temperature, T ps 1000K (e.g., Umebayashi 
1983). We also note that the realistic saturation level of the magnetic field strength is still 
unknown and remains to be determined theoretically. Obviously, this requires more detailed 
understanding of the saturation mechanism of the MRI-driven turbulence, which includes 
the understanding of the physics of rapid magnetic reconnection in the case of a high-/? 
environment (e.g., Sano et al. 2004). 

Electron Mobility Let us assume that the thermal ionization is not promising, and con- 
sider another mechanism. In weakly ionized plasma, electrons have a considerable bulk 
velocity to maintain current and, hence, magnetic field. For sufficiently high density, the 
collision frequency of electrons is larger than the gyration frequency so that the conductiv- 
ity tensor is isotropic and diagonal. Thus, we can assume that the average velocity of the 
electrons is essentially in the direction of the electric field (— v e \\E) in the comoving frame 
of the fluid. The averaged velocities of charged particles are related to the current density j 
via j = Yjitqifii (vi) or eg^ (i>j) = fcj, where the angle brackets mean a phase-space average, 
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and fi is the fraction of the current due to the ith species of charged particle. From one of 
Maxwell's equations, 4nj = cV x B, we estimate j ~ cB/(2X), where A is the typical length 
scale of the magnetic field. Thus, we can estimate the average velocities of the charged 
species: 

where we may assume f e xs 1 for (/n > 3 • 10~ 28 sec -1 cm 3 and A < 0.3H accords with the 
result of numerical simulations in Sano et al. (2004). This average electron velocity (v e ) is 
surprisingly large. Next we show its implication for the distribution function of electrons. 



Energetic Electrons To clarify the possible ionization process, we should go into the 
details of microphysics. The electron distribution function f{p) in weakly ionized plasma 
was already studied in detail (Druyvesteyn & Penning 1940; Landau & Lifshitz 1993). 

— * 

If we choose the z-direction anti-parallel to E, the electron distribution function is given 
by the following form in spherical coordinates (p, 9, 0) for the phase space: 



/*/o(p) + /i(p)cos(0), (6) 
3e 2 

fr, 

k B T-f V ~M 



/ (p) = Aexp(--^-), (7) 



AW—^/o- (8) 



where 

p 2 



2m P 



(9) 



eEl [M , 
7=— -W— »1, 10 
k B T V m e 

and m e is the electron mass and M m 3670m e is the hydrogen molecule mass. The nomaliza- 
tion factor A is propotional to the number density of electrons. From these, we can calculate 
the average energy and the average velocity in the direction of the electric field. 



W = oW?< (ii) 




(e) = 0.53M (v z y . (12) 

At first glance, it is surprising that the mass of the collision partner M appears in equation 
(12) instead of the electron mass m e . This corresponds to the very low mobility of electrons 
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in the collisional medium, which is due to m e <C M. If we insert the expression for (v e ) ~ 
42km/ sec in equation (5) into (v z ) in equation (12), we have (e) = 19eV, which is larger 
than the ionization potential of the hydrogen molecule. Thus, most of the free electrons have 
sufficient energies for ionization for this choice of parameters. This means that a significant 
fraction of the collisions may result in ionization, corresponding to a too high ionization rate. 
In reality, inelastic collisions with ionization/excitation losses should modify the distribution 
function (eq. [6]) that is obtained by considering only elastic scattering. A higher ionization 
degree results in a lower average velocity of electrons (see eq. [5]), which in turn reduces 
the ionizing collision according to equation (12). Therefore we expect that, in reality, the 
ionization degree might be on the order of x e 1CT 13 . 

5. Yet Another Mechanism to Maintain Ionization Degree in Dust-Free 

Regions 

The analysis of the ionization degree in the previous section is focused on the mid-plane 
of the disk, where external ionizing radiation is not available owing to the high column 
density of the disk. However, the gas density of the disk is a decreasing function of height 
(distance from the mid-plane), and thus, the surface layer of the disk is supposed to be well 
ionized by the incident interstellar radiation field, or by the radiation from the other stars in 
the cluster environment (see, e.g., Adams & Myers 2001). In the fully developed turbulence 
driven by MRI in the disk, the fluid elements travel from place to place according to eddies. 
The inverse cascade of turbulent energy tends to result in the state in which most of the 
power is in the largest eddies (Sano & Inutsuka 2001; Sano et al. 2004). This implies 
that each fluid element travels from the surface layer of the disk to the mid-plane, and vice 
versa. We should take into account this motion in the analysis of the ionization degree. 
For simplicity, let us consider the region with a very small amount of dust grains or the 
case in which the dust grains have grown up to larger sizes so that the total surface area 
of grains has become negligible. If we can ignore the recombination on the dust grains, the 
recombination should happen in the gas phase. The rate equation for the fractional electron 
number (x e = n e /n) is written as dx e /dt = ( — ^ . f3jX e rij = ( — /3'nxl, where f3j represents 
the coefficient for the recombination of electrons with the jth species of positively charged 
particles and (3' = ■ f3j(xj/x e ) is the effective recombination coefficient. Let us consider the 
fluid element that travels from the surface to the midplane. The approximate time evolution 
of the ionization degree x e (t) can be obtained by solving the rate equation with ( = and 
the approximately constant j3'\ 

Xe {ty l = x-l + io 11 ( j — -) ( ™ - ) ( —) , (13) 

w e '° \3- 10" 12 sec- x cm 3 J V10 15 cm ~ 3 J \lyr J V ; 
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where x e fl is the initial ionization degree. For the high-density regime without the effect 
of grains, the recombination with positively charged metals is the most important process 
(Sano et al. 2000; Fromang et al. 2002) . That is, within an eddy turnover time (~ 1 
yr at 1 AU), the ionization degree keeps the level of x e ^> 10~ 13 even without ionization 
due to cosmic rays, if we adopt the standard value of f3' ~ 3 • 10~ 12 cm 3 sec -1 for metals and 
x e fl ^> 10~ 13 . In other words, the ionized region penetrates into the neutral region by the 
turbulent motions. In this way, the turbulent motions efficiently homogenize the ionization 
degree if the recombination rate is sufficiently low. This spreading out of the weakly ionized 
region should be an important mechanism for activating MRI in the region that has initially 
low ionization. Note that this mechanism is not so effective in the region with the standard 
interstellar dust grains, because the recombination on the grain surface is much faster than 
in the gas phase. 

6. Conclusions 

Once the ionization degree becomes sufficiently high (eq.[l]), MRI-driven turbulence 
is the most promising mechanism for gas accretion in protoplanetary disks. Conversely, 
once the MRI-driven turbulence becomes active, it can consistently maintain the ionization 
degree without help of other external ionizing sources: the macroscopic energy dissipation 
in turbulence provides sufficient energy for the ionization that is required for activating 
MRI, and microscopically the energetic electrons in weakly ionized gas should provide a 
sufficient rate of collisional ionization. The turbulent spreading of the ionized gas also 
provides the mechanism of spreading the magnetorotationally unstable region. Altogether 
we conclude that most of the region in the protoplanetary disk should be sufficiently ionized, 
thus, effectively removing the dead zone. The absence of the dead zone in the protoplanetary 
disks has tremendous impact on the theoretical modeling of planet formation. 
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Fig. 1. — Fractional number densities of electrons (e _ ), metals (M + ), and molecules (mo + ) 
as a function of the ionization rate divided by the total number density in the model with 
dust grains. Recombinations on grain surfaces are taken into account as well as the radiative 
and dissociative recombinations in the gas phase. The size of the dust grain is assumed to 
be 0.1 (Km. See Sano et al. (2000) for details. 
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Fig. 2. — Schematic diagram for the feedback loop in the turbulence driven by magnetorota- 
tional instability (MRI). Once the protoplanetary disk is sufficiently ionized, the MRI occurs 
and the magnetohydrodynamic turbulence develops and saturates. At the saturated state, 
most of the energy dissipation results in the heating of the gas and its conversion into the 
radiation that eventually escapes from the disk. If only a small fraction of the energy dis- 
sipation is used for ionization, the disk keeps this magnetic activity (see eq.[4]). One of the 
possible processes of ionization should be the collision between the energetic electrons and 
neutral particles (see Section 4). 



